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I. INTRODUCTION 



In recent times major advances in the experimental techniques have led to the realization of experiments in which 
quantum systems in a single particle regime are studied under their controlled interaction with some environment. A 
paradigmatic example in this context is given by the motion of a massive test particle in an interferometric setup, which 
gives rise to interference fringes as typical quantum signatures. When the coupling with the environment becomes of 
relevance, such interference fringes are gradually washed out, and a classical dynamics is eventually recovered. This 
phenomenon goes under the name of decoherence [TH3]. Its understanding and theoretical description require on the 
one hand a control over the environment, on the other hand a microscopic model for the interaction and the ensuing 
dynamics. 

For the case of a tracer particle immersed in a dilute gas such a microscopic description has been obtained considering 
the centre of mass degrees of freedom only. The reduced dynamics is given by a master equation in Lindblad form 
which has been called quantum linear Boltzmann equation, since it provides the natural quantum counterpart of the 
classical linear Boltzmann equation (see [4j for a recent review and references therein). The microscopic input is 
given by the complex scattering amplitudes describing the collisions between gas and test particle, while the gas is 
characterized by its density and momentum distribution. In this paper we consider an extension of this result, which 
includes internal degrees of freedom of the tracer particle. The microscopic derivation is performed along the lines 
of a general strategy for the derivation of Markovian master equations, which relies on a scattering description of 
the interaction events |5j. Besides the gas properties, this approach takes as basic input the multichannel complex 
scattering amplitudes, which describe the influence of the internal states on the scattering events. Indeed, when 
the scattering cross section does not only depend on the relative motional state between tracer and gas particle, 
such an extension becomes mandatory in order to correctly describe the dynamics. According to the Markovian 
approximation, the obtained master equation is in Lindblad form. This derivation confirms the structure of the 
dissipative term, which has been heuristically obtained in [B], further determining the coherent contribution to the 
dynamics due to forward scattering. The latter becomes relevant in the determination of the index of refraction for 
matter waves. When either type of degrees of freedom can be described in classical terms, a Markovian quantum 
classical master equation is obtained. Such a result, corresponding to a classical treatment of the motional degrees 
of freedom, has been considered in [7]. In that context the name Bloch-Boltzmann equation was proposed for the 
equation, since for a two-level system an extension of the optical Bloch equations to include a Boltzmann-like collision 
term is obtained. In the same spirit, the name quantum Bloch-Boltzmann equation can be used to indicate a master 
equation, which gives a quantum description of both internal and centre of mass state. 

An interesting situation appears when in the final detection the internal state of the test particle is not resolved 
at the output of the interferometer. In this case the internal degrees of freedom become part of the environment. 
Then a non-Markovian dynamics for the motional state appears, which can be described in terms of a coupled set of 
Lindblad equations for the unnormalized statistical operators corresponding to specific internal channels. This type 
of non-Markovian dynamics can be considered as a generalized non-Markovian Lindblad structure. It arises as a mean 
over a classical index, which can take place e.g. as a consequence of the interaction with a structured reservoir [6l|8l[9]. 
This situation is here considered in the study of the loss of visibility of the interference fringes in an interferometric 
setup. The ensuing decoherence effect is generally not described as an exponential loss of visibility depending on the 
strength of the interaction, as in the usual Markovian case. 

The paper is organized as follows. In Sect. |TT] we consider the expression of the master equation, pointing to the 
main steps necessary for its derivation and putting into evidence the microscopic quantities determining its explicit 
form. A detailed microscopic derivation of the master equation is performed in Appendix |A] The master equation is 
given both in terms of matrix elements of the statistical operator in the momentum and internal energy eigenstates 
basis, as well as an explicit operator expression, which makes its Lindblad structure manifest. This also allows to 
easily recover under suitable limits previously considered master equations, which describe either only one of the 



two kind of degrees of freedom or a hybrid quantum classical description of both. In Sect. Ill we show how the 
interplay between internal and motional states can influence the visibility in an interferometric setup for the study of 
decoherence, leading to a non-Markovian behaviour in the reduction of the visibility of the interference fringes. 



II. THE MASTER EQUATION FOR A TEST PARTICLE WITH INTERNAL DEGREES OF FREEDOM 

We first consider the key ingredients and steps which lead to obtain the master equation describing the coUisional 
dynamics of a test particle immersed in a structureless background gas, keeping the internal degrees of freedom of 
the particle into account. The task of a full microscopic derivation will be accomplished in Appendix [Xj relying on 
a method recently introduced for the derivation of Markovian master equations, which has been called monitoring 
approach [SJ llOH12j . In the monitoring approach the reduced dynamics of a system in contact with some environment 
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is obtained describing their interaction by means of scattering theory. The building blocks in such a formulation of the 
open system dynamics are therefore the S-matrix characterizing the single interaction events and the rate of collisions. 
Both quantities are given by operators on the tensor product Hilbert space of system and environment, which we shall 
denote by S = I + iT and f respectively. The operator nature of these quantities is crucial in order to keep the gas 
and test particle state into account in the dynamic description of the coUisional interaction. The Markovian master 
equation for the reduced dynamics is obtained by assuming the various collisions as independent, so that their effect 
cumulates according to the state dependent scattering rate, and taking the trace over the environmental degrees of 
freedom. 



A. Expression of the master equation 



The formal expression of the master equation reads [S] 



(1) 



where H is the free Hamiltonian and p is the statistical operator of the system. For the case at hand the free 
Hamiltonian of the system is given by 



H = 



2M 



(2) 



where P is the momentum operator of the test particle, M its mass and n basis of energy eigenstates 

in (D". The superoperators appearing at r.h.s. of Eq. ([T]) are defined according to 



and 



7^p = ^Trgasf^^l/2 Re (T) T^/^ 



(3) 



(4) 



respectively, where Pgas is the single particle statistical operator describing the gas environment. Note that the 
operators C and TZ arise by acting with an operator in Lindblad form on a state of system plus gas in factorized 
form, further taking the partial trace with respect to the gas. While this operation is formally legitimate, and 
guarantees preservation of trace and Hermiticity of the statistical operator describing the test particle, it is generally 
not true that the resulting dynamics for the reduced system only is given by a master equation in Lindblad form, 
thus granting complete positivity and describing a well-defined Markovian dynamics. Indeed this step involves further 
approximations, which depend in a crucial way on details of system and interaction. It is well known that by taking 
the partial trace with respect to the unitary evolution of the overall system one can obtain a Markovian dynamics only 
if further hypotheses hold. This remains true for the case at hand, despite the fact that important approximations 
have already been introduced in replacing the Hamiltonian dynamics for system plus gas with a Lindblad operator 
only specified by T and f. The actual proof that a Markovian dynamics applies to the situation of interest and 
the specific expression of the superoperators appearing in Eq. ([T]) is obtained through the microscopic calculations 
performed in Appendix \K\ 

Relying on the results of Appendix Ala we write the following expression for the contributions in Eq. ([3| in the 
momentum and channel basis i)} 

{P,i\Cp\P',k)^Yl J '^Q [{P-Q,MP' -Q.l)Mg{P,P';Q) 



--{P,j\p\P', k)Mii{P + Q,P + Q-Q) 
- i (P, z|p|P', l)Mf^{P' + Q,P' + Q; Q) 



(5) 



4 



where the complex rate functions Mi^ {P, P'; Q) are given by 

Mg (P, P'- Q) = x^I / dp {p, P-Q;Q) Ll^ {p, P' ~Q;Q), (6) 

with the Xik ^ notational shorthand to indicate that the contribution is different from zero only for £ij — £ki , where 
Skj = Ek — Ej denotes the difference in energy between internal states, while the p-integration is restricted to the 
plane Qj^ = {p e R'^ : p • Q = 0}. The functions Lij {p, P; Q) are defined according to 



x/i, ( rel(pj.,P±) - § + 7i^Q.t<HpL,PL) + § + T^irQ 



(7) 



where denotes the stationary gas distribution, rtgas is the density of the gas, m the mass of the gas particles, 
m* = mM/{M + m) the reduced mass, and fkj{Pf,Pi) denote the multichannel complex scattering amplitudes, which 
depend on the microscopic interaction potential and describe scattering from an initial momentum pi and internal 
state j to a final state with momentum pf and internal state k. Moreover, P± and indicate, respectively, the 
perpendicular and the parallel component of the momentum P with respect to the vector Q; while rel(p, P) = 
{m^,/m)p — {m^,/M) P is the the relative momentum between the gas particle momentum p and the test particle 
momentum P. 

Exploiting these results we can easily write the master equation Eq. ([I]) directly in operator form. In fact using the 
functions (p, P; Q) let us introduce the following family of jump operators 

Lq,p,£ = e»« >^/'"' J2 ^'^0 (P> P; Q) ® E„-. (8) 

where X and P are position and momentum operators of the test particle, while the operators E^j = act on the 
internal degrees of freedom only, since |z) denotes the energy eigenstate with eigenvalue huii, and the exponential factor 
describes momentum exchanges according to exp (iX • Q/h)\P) — \P + Q). Note that the functions {p,P;Q) 
essentially depend on the scattering amplitudes and the momentum distribution of the gas, thus keeping into account 
all the details of the collisional interaction. These expressions appear operator-valued in the master equation, being 
evaluated for P — > P, so as to take into account the actual momentum of the colliding test particle. The incoherent 
contribution C in Eq. ([T]) finally reads 



(9) 



The superoperator TZ of Eq. Q on its turn according to Appendix Alb amounts to the commutator with an effective 
Hamiltonian given by 

H„ = ~2nh^'^ V / dpoAi(po)Re [/y (rel (po, P) , rel (po, P))] ® E,,. (10) 



Using the alternative expression given by Eq. ( A18[ ) to define Lij (p, P; Q), it is immediately seen that the incoherent 
term of this master equation confirms the result tieuristically obtained in [6. . In the latter reference this equation 
has also been termed quantum Bloch-Boltzmann equation in that it provides a quantum description of both motional 
and internal degrees of freedom, thus extending the result of |13j , where the centre of mass degrees of freedom were 
treated classically and the name Bloch-Boltzmann equation was used. These names should not confuse the reader. 
Indeed, only for the case of an atom in a two-level approximation undergoing a collisional dynamics this equation 
refers to an extension of the optical Bloch equations with a Boltzmann collision term. 



B. Limiting forms 



As a compatibility check of the master equation derived in Appendix [Xj and in order to make contact with previous 
work, we will now show how in suitable limits it recovers already known equations. Since the equation describes the 
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quantum dynamics of a test particle with both internal and translational degrees of freedom immersed in a dilute gas, 
natural limiting situations appear considering a structureless test particle or an immobile system. These situations 
correspond to the quantum linear Boltzmann equation |12j . and to the master equation for an immobile system 
interacting through collisions with a background gas [51 114j . Another natural limit consists in a hybrid quantum 
classical description, in which either internal or centre of mass degrees of freedom are treated classically. The master 
equation corresponding to this last case has already been considered in [13]. A classical treatment of both kinds of 
degrees of freedom leads to the master equation for a classical Markov process, with a probability density depending 
on both a discrete and a continuous index. 

a. Quantum linear Boltzmann equation If the internal degrees of freedom can be disregarded the sum in Eq. ^ 
has a single non vanishing contribution, so that instead of the multichannel scattering amplitudes fij(jJf,Pi) there is 
a single amplitude which can be indicated as f{Pf,Pi)- The incoherent term in the master equation reduces to 



(11) 



with 



X / (^rel (p^ , ) - I , rcl , Pi ) + I 
while the Hamiltonian term reads Hq + H„, where Hq — P^/2M and 



H„ = ^ / dpofi{po)Re [f (rel (po, P) , rel (po, P))] 

J 



takes into account the energy shift due to forward scattering. This result complies with the quantum linear Boltzmann 
equation obtained in [12], whose properties have been discussed in detail in 

b. Immobile tracer particle We now consider the opposite situation, corresponding to an infinitely massive test 
particle, so that the dynamics of the translational degrees of freedom can be neglected. To consider this limit it is 
convenient to come back to the expression Eq. ([5]) of the quantum master equation in terms of the complex rate 
functions mH (P, P';Q), which in the limit M — > oo when integrated over Q reduce to 

' - —)di I dp I dpoS E,>j /.,(p,Po)A!;(p,Po), (12) 

where no dependence on the test particle's momentum is left. The matrix elements of the incoherent part of the 
quantum Bloch-Boltzmann equation are therefore now given by 

mcp\k) = J2 ({j\p\i)M^i ~ \mmii - \m)M^\] , (13) 

while the coherent part corresponds to a effective Hamiltonian whose matrix elements in the energy eigenbasis are 
given by E]i = —2Trfi^{ng^s/m^)xi^ J c?PoM(Po)R-e [fij (po,Po)], thus confirming the result obtained in [5^ for the case 
of a non degenerate Hamiltonian. 

c. Quantum classical description The limiting expressions of the quantum Bloch-Boltzmann equation, obtained 
when either the internal or the translational degrees of freedom can be treated as a classical label, correspond to 
hybrid quantum classical descriptions, which naturally arise when decoherence affects on different time scales the two 
kind of degrees of freedom. 

When the centre of mass degrees of freedom can be treated classically it is convenient to introduce the classical 
rates 

M^l{P + Q;Q):=M^l{P + Q,P + Q;Q) (14) 
with (P, P'; Q) as in Eq. so that the semiclassical Bloch-Boltzmann equation reads 



d /T,N 1 



^;i(^,K)(z| + H„(P),p(P) 

ijkl 

-^Mr^(P-KQ;Q){ELE,„p(P)} 



J2 J dQ[M^(P;Q)E,,p(P-g)E 



(15) 
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where H„(P) is obtained from Eq. ( |10[ ) with the replacement P — > P, and p{P) denotes a collection of trace class 
operators in (D" normalized according to / dPTrc;™ p{P) — 1- 

If the classical approximation applies for the internal degrees of freedom the incoherent term of the master equation 
giving a quantum description of the translational dynamics only takes the form 



dt' 



1 

ih 



2M 



H! 



dp 



(p, P; Q) p,U, (p, P; Q)^ e^^^' 



|{L,, (p,P;g)^iy (p,P;Q),p,} 



Here 



H! 



2 "gas 
TO* 



c'Po/^(Po)R.e (rel (po, P) , rel (po, P))] , 



(16) 



(17) 



and Pi denotes a collection of trace class operators in L^(R"^) normalized according to X]"=i "^"^L^i^^) Pi ^ 1- 

For the case in which all the off-diagonal elements with respect to momentum and internal energy eigenvalues 
vanish, that is to say (P, i|p|P', fc) = if P 7^ P' or i fc, the motional state of the test particle is fully characterized 
by the distribution of the diagonal terms fi{P) = (P, z|p|P, i), which is a classical probability density obeying the 
classical Markovian master equation 



^/i(P) = ^ 1 dQf,{P - Q)Mli{P- Q) - /,(P) ^ j dQMf^iP + Q- Q), 



(18) 



where the positive quantities Mf^ (P; Q) defined in Eq. ( 14 ) can actually be interpreted as the transition rates from an 
initial momentum P — Q and internal state j to a final momentum P and internal state i. This classical Markovian 
master equation provides the natural generalization of the classical linear Boltzmann equation to a particle with 
internal degrees of freedom |15| . 



III. EFFECT OF INTERNAL DEGREES OF FREEDOM ON CENTRE OF MASS DECOHERENCE 



The quantum linear Boltzmann equation has proven useful in the description of collisional decoherence, as well 
as in the evaluation of the index of refraction for matter waves [H [T^l [TBHUj • We will now consider the effect of 
internal degrees of freedom, affecting the collisional interaction between massive test particle and background gas, on 
the visibility of the interference fringes in a interferometric setup. In particular we will show that the visibility can 
exhibit oscillations due to non-Markovian effects. The effect of the entanglement between internal and centre of mass 
degrees of freedom for the visibility of quantum interference experiments has already been considered in |19) , in the 
absence however of decoherence effects. 



A. Generalized Lindblad structure 



The quantum master equation Eq. ([T]) is in Lindblad form: this means that the dynamics of the test particle is 
Markovian when both translational and internal degrees of freedom are described and detected. A different situation 
emerges if the translational or the internal degrees of freedom, although influencing the collisional dynamics, are not 
revealed during the measurement process. In this case they must be averaged out from the description of the system, 
by means of the partial trace, thus becoming part of the environment. As well known in the classical case, a non- 
Markovian dynamical regime becomes Markovian by suitably enlarging the set of degrees of freedom and vice-versa. 
Indeed, a unitary Markovian time evolution for both system and reservoir generally gives a non-Markovian reduced 
dynamics for the system, the degree of non-Markovianity of the description also depending on where we set the border 
between system and environment, which ultimately depends on the physical quantities actually measurable by the 
experimenter. A smaller set of observed degrees of freedom, with respect to those actually involved in the dynamics, 
can lead from a Markovian to a non-Markovian regime. A general mechanism describing this passage in quantum 
systems is presented in : a Lindblad structure on a bipartite system can generate in the two reduced subsystems a 
generalized Lindblad structure, typically describing a non-Markovian dynamics. 

In the situation we are considering, the bipartite system is formed by the translational and the internal degrees of 
freedom of the test particle. If the measurements at the output of the detector cannot probe the internal degrees of 
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freedom, the only experimentally accessible quantities are expectations or matrix elements of the statistical operator 
given by 



2it) ^ Trc. {pit)} = ESPm) E/^^W' 



(19) 



where p{t) is the statistical operator describing the full dynamics of the test particle. It is easy to see that, if the free 
Hamiltonian is non degenerate, the diagonal matrix elements in the energy basis with respect to the internal degrees 
of freedom of the master equation lead to Eq. (16), i.e. a coupled system of equations for the collection {pi{t)}- of 
trace class operators on L^(R'^). This system oiKiuations has a generalized Lindblad structure [5], and therefore it 
can also describe highly non-Markovian dynamics for the statistical operator g{t) given by Eq. ( |19[ ). Indeed, there is 
generally no closed evolution equation for g{t), but from the knowledge of the initial collection {p i(0) },- the generalized 
Lindblad structure allows to obtain the collection {pi{t)}^ at time t and therefore, through Eq. (19), also g{t). In the 



next paragraph we are going to explicitly point out non-Markovian behaviour described by the generalized Lindblad 
structure, which express the effect of correlations between internal and translational degrees of freedom on the visibility 
of interference fringes for superpositions of motional states. A complementary situation has been considered in [20] . 
where the effect of coUisional decoherence on internal state superpositions of a cold gas has been studied in detail. 

In typical interferometric experiments the test particle is much more massive than the particles of the background 
gas. The dependence on the momentum operator in the Lindblad operators describing the collisional dynamics and in 
the Hamiltonian part determining the energy shift can therefore be replaced by a fixed value Pq, which represents the 
initial momentum of the test particle entering the interferometer. Taking the diagonal matrix elements of the general 
form of the master equation given by Eq. ([I]) and specified by Eq. (10) and Eq. (|9]), assuming non degeneracy of the 
internal energy eigenvalues one finally obtains for the collection of operators Pi(t) = {i\p{t)\i) the following coupled 
system of equations: 



dt 



P.it) 



1 

ih 



2M 



,Ptit) 



iQX/h 



it) 



(20) 



with Vp and Tp probability densities and transition rates defined by 



iQ) ■= 



MlliP^Q) 



JdQMl^iP + Q;Q) 



and 



T% :^ J dQM^iiP + Q;Q). 



(21) 



(22) 



respectively. Note that at variance with Eq. ( 16 1 we are now not assuming a classical dynamics for the internal 



degrees of freedom, we focus on the diagonal matrix elements of the internal states only since the latter are enough 
to determine the non-Markovian dynamics of the motional state according to Eq. (19). The fact that the positive 



quantities (P; Q) are transition rates implies that Vp iQ) can be interpreted as the probability distribution 
function for a test particle with momentum P and internal energy eigenstate j to exchange a momentum Q, and 
to go into the internal state i due to a collision with the gas. On the same footing, Tp can be interpreted as the 
total transition rate for a test particle with momentum P and internal state j to go to a fixed final internal energy 
eigenstate i. 



B. Explicit solutions in position representation 



We are now going to describe the visibility reduction predicted by the generalized Lindblad structure Eq. (20) 
obtained from the quantum Bloch-Boltzmann equation Eq. ([I]) in the limit of a very massive test particle. 

To obtain the formula describing the fringes visibility in an explicit way we need to solve the equation of motion in 



the position representation. Starting from Eq. ( 20 ) and omitting for simplicity the explicit dependence on the classical 



label Pq denoting the momentum of the test particle, we obtain 

-^p,(X, X', t) = ^(Ax - ^X')p^iX, X', t) 
at in 



E (r''*'' iX - X') (X, X',t)- PV.(X, X', t)) , (23) 
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where pi{X , X' ,t) denotes the matrix element {X\pi{t)\X') and {X — X') is the characteristic function of the 
probability density V^^{Q) [21], i-e. its Fourier transform 



{X -X')^ J dQe'(^-^'' '?/'"T'^ (Q). 



(24) 



We will now consider a few cases in which Eq. ( 23 1 can be solved analytically, so as to obtain an exact expression for 
the visibility, showing up different possible qualitative behaviour. 



1. N-level system 

When the collisions are purely elastic, so that they do not lead to transitions between different internal states, the 
scattering rates satisfy F*^ = iJjjF". This is the case when the energy exchanges involved in the single collisions are 
much smaller than the typical separation of the internal energy levels [13) . The equations for the different pi then 
become uncoupled and take the form 

^o^iX, X', t) = ^(Ax - ^X')P^{X, X', t) - F"(l - {X - X'))p,{X, X', t). (25) 
dt in 

The latter equation can be conveniently solved introducing the function |22| 

X.(A, M, t) := Tr {p,(t)e'(^ >^+'^ P)/''} (26) 



where X and P as usual denote position and momentum operators of the test particle. In such a way Eq. (25 1 leads to 

A 



• -r"(i -$"(/.)) 



X.(A,/x,i), (27) 



which is an equation of first order solved by 

X.(A,/x,<) =X°(A,Ai/Af + /x)e-r"^o(i-*"Wt-0/A^+M))dt'^ (28) 



where the function x?(A, Ai/M + /x) obeys the free equation 9tXi(A, /i., i) — (A/M) • 9^Xi(A, /x, t). Inverting Eq. (26) 
by taking the Fourier transform with respect to A, 

p.{X,X',t) = J ^^e-^-(^+^')/2\,(t,A,X-X'), (29) 

we obtain the exact solution 

pdX,X',t) = f ,^^e-'^-^/''e-r"/o'(i-*"(A(t-0/A/+x-x'))dt'^0(^^^^^,^^^^) 

J [2TThy 

expressed in terms of an integral of the freely evolved subcoUections p°(X,X',t) with a suitable kernel, where we 
have set 



pUX,X',t)^ J ^e-^ (^+^')/^V(A,X-X',t) 



and pi{0) — {i\p{0)\i). The evolution of the statistical operator given by Eq. (19) is obtained summing the different 
Pi{X, X' ^ t) over the discrete index i. For an initial state given by a product state between the translational and the 
internal part, so that Pi{Q) — PiQiO), we finally obtain 

g{X, X', t)^ypj -^^e-^^-^/^e-i^" /„'(i-$-(A(t-t')/M+x-x')) dt'^o^;^ + s,X' + s, t). 

^ J (2'Kny 

(31) 

This result reduces to the standard Markovian situation, when either only one Pi is different from zero (and therefore 
equal to one), or the rates are all equal. This limiting cases describes situations in which the initial state is in a 
specific internal state or the collisions do not depend on the internal state of the tracer particle. 
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2. Two-level system 



For the case of a two-level system a natural situation corresponds to inelastic scattering taking place only when the 
test particle gets de-excited, so that only one of the two scattering rates is different from zero. This case can still be 
treated analytically. Assuming F^^ = 0, the equation for X2{t,\,fi) gets closed, and is solved by 

X2(A, M, t) = x^(A, Xt/M + M)e-r"te-r" /,'(i-$-(A(t-t')/M+M)) dt'^ 

The equation for xi(A, then reads 

9tXi(A,/i,t) = ^ 

and its solution is given by 

Xl(A, t) = e-^" /,:(l-*^^(A(t-*')/M+^))dt' |^0(;^^ ^^/^j ^ ^) 



9^-F"(l-a>"(/x)) 



Xi(A,/x,i)+F^^$^^(/x)x2(A,/x,i) 



(32) 
(33) 



+T 
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,r"/n il-<S>'\X{t^t")/M+fj.))dt" 



'.<^>^^{\{t - t')/M + At)x2(t', A, \{t - t')/M + ii)]dt' I . 



(34) 



This formula explicitly shows that xi(A, /.t,i) depends on the function X2(A,/x, •) evaluated over the whole time 
interval between and t, a typical signature of non-Markovian dynamics. Assuming once again that the initial state 
is characterized by Pi(0) — Pig{0), the statistical operator describing the translational degrees of freedom of the test 
particle is given at time t by the expression 



P2e ■'^"tf.-T^^^ Io(^-'^^^W^-t')/^''+^'^'))'i^' 



+p^e-^'' foi^-'i'''(Mt-t')/M+x-x'))dt' ^p_^pi2g-r"/„*(i-<t."(A(t-t')/M+x-x'))dt' 

|^g-r^'t'g-r^Vo\i-*''(M«-t'')/J>/+x-^'))dt''gr"/j'(i-*"(A(t-t'')/M-HX-x'))dt'' 



x$i2(A(i - t')/M + X -X') 



dt' 



(35) 



C. Non-exponential visibility reduction 



We can now explicitly present the visibility reduction predicted by the generalized Lindblad structure obtained 
from the quantum master equation for a test particle with internal degrees of freedom. Our aim is to obtain an 
exact expression for the loss of visibility in a double-slit arrangement as a function of the time of interaction with the 
environment, and to illustrate by means of example how the presence of the various scattering channels, corresponding 
to the different internal states, can actually lead to non-Markovian behaviours. In particular we will consider the 
situation of purely elastic collisions in full generality, also allowing for inelastic scattering in the case of a two-level 
system. While the experimental setting is always taken to be the same, the different number of internal degrees of 
freedom involved and the presence or absence of inelastic scattering events will lead to more or less marked non- 
exponential behaviours in the reduction of the visibility fringes. 



1. Visibility formula 

We first derive a formula for the visibility reduction in the case of a double-slit experiment in the far field approx- 
imation. A beam of particles moves towards a grating perpendicular to its direction of propagation, and with two 
identical slits separated by a distance d, finally reaching a detector where the fringes of interference are observed. 
During the flight through the interferometer the beam particles interact through collisions with the environment in 



the background, thus undergoing decoherence. We consider an initial product state, so that in the notation of Eq. ( 20 ) 
one has pi(0) — Pig{Q), where g{0) describes the translational degrees of freedom. If after the passage through the 
collimation slits the test particle is described by psi, then the double-slit grating prepares the initial state 

e(0) — 2 cos psi cos . (36) 



2h / \ 2h 
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Setting 



we consider the quantity 



giX,X,t)^{X\Uit)[gm\X) := I (X) 



V = 



which describes the reduction of the interference pattern with respect to the free case. Exploiting the fact that the 
time evolution generated by Eq. (31 1 is covariant under translations [4l[24], so that 



mi' 



one has, using Eq. (361, 



2|(X-id|^^ft)[p3ie-'P-^/^]|X-id)| 
[X - ld\Uit)[p,i]\X - \d) + (X + \d\U{t)[p,i]\X + \d) ' 



(37) 



(38) 



where now t is the time employed by the test particle to reach the detector. Indeed this result remains true for any 
translation-covariant time evolution. 



For an initial factorized state of the test particle we can exploit Eq. (31) to obtain a closed formula for the time 



evolution operator lA{t) depending on the initial internal state, i.e. on the coefficients Pi appearing in Pi(0) = PiQ{0): 
the numerator of Eq. ( 38 ) then reads 



dsdX ^ — iX- 



'/^{X - f + s|Wo(t)[psie-^P"'/'']|X- f + s)e-r'Vo(i-*"(^(*-*')A^))dt'|^ 



where U^ii) is the free evolution operator of the translational degrees of freedom, so that 



(X- 



d 



s|Z^o(t)[Psie-*P"^/'']|X- 



The latter expression can also be written 
(X 



^|z^o(t)M|x + ^) 



^) = (^-f + ^|Wo(t)M|x + ^ 



(39) 



xe 



^iMd-X/(ht) / dYdY iM(Y^-Y''^)/(2ht) 

J [2nhf 

iMX-{Y-Y')/{ht)^iMd-{Y+Y')/(2ht)^Y\p^^\Y'), 



(40) 



assuming due to symmetry Tr (Xpsi) — 0. 

This formula enables us to implement the far field approximation. In fact, let cr be the width of the two slits, so 
that the integrand is negligible if Y (and similarly for Y') takes values outside the support of ps\, then MY"^ / [ht) < 
Ma^ I (ht) and therefore for a time long enough such that ht/M ^ cr^ the first exponential can be disregarded. The 
same applies for the last exponential if ht/M ^ ad. For times longer than max^AIa^ /h, Mad/h^ , corresponding to 
the far field approximation, we get 

{X imt)[ps^]\X + ^)^( ^Xe-^Md.x/m~^^ ( M\ ^ ^^^^ 



where psi (•) is the distribution function for the momentum of the particle in the state ps 



Psi I —X 



dYdY' 



-.MX-iY-Y')/iM)i^Y\p,i\Y' 



The equivalence between the assumption ht/M ^ a"^ and the far field approximation i ^ cr^/A, where A — h/P^ is 
the wavelength associated to the test particle and L is the distance between grating and detector, is easily seen from the 
relation L = p^t/M , where Pz is the component along the z direction of the massive particle, assumed to be constant. 
Substituting Eq. (39 1 in the numerator of Eq. (38 1 and using the approximation psi {M{X + s)/t) w psi {MX /t) valid 



because of the localization of the state ps\, we can easily perform the integrals over s and A, thus finally obtaining 



t 



For the denominator of Eq. ( |38[ ) one can proceed in an analogous way, using 

' M^ 3 
T 



(X ± + s|Z^o(i)M|X ± + s) 



Ps\ 



M 



and performing the integral over A, further observing that (0) — 1 for the normalization of V'^^{Q). 



11 



2. N on- exponential behaviours 



The desired expression for the visibiHty in the absence of inelastic scattering and for an arbitrary number n of 
channels thus reads 



V 



(42) 



where we recall that the probabilities Pi give the weight of the different internal states in the initial preparation. The 
dependence on t in this formula can be easily made explicit with the change of variable t' /t = s, so that one has 



(43) 




FIG. 1. Plot of the visibility in a double-slit arrangemen t as a function of the interaction time with the environment, for 
the case of elastic scattering events only, according to Eq. (431 and with growing number of channels from left to right. The 



dashed lines represent the Markovian exponential decays occurring if a single elastic channel prevails on the others, the one 
with the highest and lowest decay rate corresponding to lower and upper line respectively, (a) Visibility for n — 2 elastic 
channels, according to the expression Eq. (461. It appears a non monotonic decay as a consequence of t he interference between 
the contributions of the two different elastic channels. The coefHcients a' and /3' defined in Eq. (45 I are calculated for two 



Gaussian distributions V (Q) and V (Q) of the exchanged momenta. Taking d = dz as direction of propagation inside the 
interferometer we only need to specify the mean and the variance of the exchanged momenta along this axis, respectively j-iu 
and an, i = 1,2. The plot is for pi = P2 = |, while F^^ = F^'^ = 10, d = 1, an = (T22 = 0.1, fin — —0.2, fi22 ~ 0.3 in arbitrary 
units, (b) Visibility for n = 8 elastic channels according to the general expression Eq. (43 1. The characteristic functions "1>" 
are calculated starting from Gaussian distributions, assuming equal rates F" — 10 and equal variances an =0.1 in arbitrary 
units as in (a). The pi are uniformly distributed and the means fiu are equally spaced in the range from —0.2 to 0.3. 



From Eq. (43 1 one can easily see the difference between the Markovian situation, corresponding to n = 1, and the 



general case. If there is just one ter m in the sum, the modulus simply picks out the real part of the characteristic func- 
tion in the exponential and Eq. (43 ) describes an exponential decay in time with a rate F ^1 — Re {$ {d{s — 1))} ds^ . 

This can happen if only one internal energy state is populated in the initial preparation or the scattering events are 
actually independent on the internal state. If there are at least two terms, the modulus can generate oscillating terms 
as a consequence of the interference of the different phases arising since the functions are in general complex 
valued. Even if the imaginary parts of the characteristic functions are zero, i.e. the distribution functions of the 
exchanged momenta are even, Eq. (43) can describe highly non-exponential behaviour. In this case in fact it reduces 
to 



i=l 



'{l-Ii<i'"{d(s-l))ds)t 



(44) 



i.e. the sum of different exponential functions. As shown in [51 [5S] this kind of relations can describe behaviour very 
different from the exponential one. 
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Let us consider in more detail the case of a two-level system. Introducing the notation 



a' :=Rc / {$**(d(s- l))}ds 

:=Im f {$"(d(s-l))}ds, 
Jo 



(45) 



the visibility reduction is explicitly given by 
V 



2 ~2V^^(l-o^)t I 2 -2r^^(l-Q^)t , o -r"(l-QMt -r^^(l-a^)t „„„ r/-pll ol -p22o2\,l 

p-^e ^ ' + ^ ' + 2pip2e ^ ' e ^ ' cos [(1 p — 1 P )t\ 



1/2 



(46) 



This formula describes a decrease modulated by the oscillations produced by the cosine function. To illustrate this 
behaviour in Fig. [T]a we plot the visibility as a function of time, considering by means of example two Gaussian 
distributions. Note that the a ppea rance of the oscillations depends on a non vanishing mean value for the distribution 
functions V^^{Q) given by Eq. ( 21 ), which describe the state dependent momentum transfers. This feature corresponds 
to a preferred direction in the net momentum transfer between test particle and environment, as happens e.g. by the 
interaction with a laser beam [26J, the asymmetry in the single interaction channel being determined in this case by 
the direction of propagation. 

The behaviour described by Eq. (43 1 for an n-level system is illustrated in Fig.jljb, where we show how the increased 
number of levels can strongly suppress the oscillations and lead to a reduction of the visibility. The dashed lines 
represent the exponential decays pertaining to the Markovian situation arising if only one of the internal energy states 
is initially populated, the one with the highest or lowest decoherence rate corresponding to the lower or upper dashed 
line respectively. It appears that with growing n the interference between the contributions of the different channels 
to Eq. ( 43 ) rapidly determines a decay of the visibility sensibly faster than that occurring for the corresponding 
Markovian single-channel dynamics. Indeed in Fig. [T] left and right panel correspond to the same interaction strength 
but differ in the number of involved degrees of freedom, ranging to n = 2 to n = 8. 

Relying on the results of Sect |IIIlT2| one can also obtain an expression of the visibility in the presence of inelastic 
scattering for a two-level system. Indeed starting from Eq. (35) and following the same procedure as above one comes 
to 



V = 



pl2g-riit/„i(l-<[.ii(d(^-l))ds 



't'e+r''/o'(i-*"(d(t"-t)/0)dt"^i2 / ^ 



,t' -t 



dt' 



(47) 



where for simplicity P2 = 1, and we have taken F^^ = 0, so that the oscillations in the visibility cannot be traced back 
to interference among different components. An illustration of the behaviour of the visibility in this case has been 
plotted in Fig. [2] always assuming for the sake of generality a Gaussian distribution of momentum transfers. In this 
case the dashed line corresponds to the exponential Markovian decay occurring if only the elastic channel is involved 
in the dynamics. It immediately appears that a non monotonic behaviour in the loss of visibility is observed also in 



this case, due to the multiple time integration in Eq. ( 35 ) 



We have here considered the visibility reduction as a function of time. However in typical interferometric experiments 
the time of flight is fixed, and it is more natural to study the loss of visibility as a function of the strength of 
the interaction with the environment. In coUisional decoherence this depends on the number of collisions, directly 
proportional to the gas density or equivalently to its pressure for a fixed time of flight. As in the Markovian case, 
we can thus express the visibility as a function of the pressure of the background gas, which is the physical quantity 
directly tunable in actual experiments [T5]. Introducing the effective cross section (7es{Po, i) according to the relation 

m 



E 

3 



dQMfAP + Q-Q) 



'M 



where the l.h.s. denotes the classical loss term appearing in Eq. (18), one has for an ideal gas 



Foo-cff(-Po,*), 



(48) 



where p is the pressure of the gas and T its temperature. One can thus introduce a family of reference pressures 
depending on the initial internal state of the particle entering in the interferometer 



MknT 



(49) 
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FIG. 2. Plot of the visibility in a double-slit arrangement as a function of the interaction time with the environment, for the 
case in which one of the internal states also undergoes inelastic scattering, according to Eq. | |47[ ) with n = 2. It clearly appears 
a non monotonic decay of the visibility as a consequence of the multiple time integration describing the contribution of the 
inelastic channel. The distributions of momentum transfers are assumed Gaussian, with tJn = l,a"i2 = 3, /in — l,/ii2 ~ 5; 
moreover F^^ — 0.75 and F^^ — 1.75. The dashed line corresponds to the Markovian dynamics determined by the channel 
undergoing elastic scattering only. 



where t is the time of flight, so that Eq. (43 1 can equivalently be written as a function of the pressure in the 
interferometer 

r^o^ = 4- (50) 

Po 

This simply imphes that the behaviour of the visibility as a function of time is equivalent to its behaviour with respect 
to the pressure and therefore the disturbance of the environment. 



IV. CONCLUSIONS 



We have derived the master equation describing the dynamics of a test particle with both translational and internal 
degrees of freedom, interacting through collisions with a low density background gas. This has been done building on 
the so-called monitoring approach [S] , and confirms a previous heuristic argument put forward by one of us in [B] . The 
present microscopic derivation further allows to determine the energy shift. As we have checked, the result reduces 
to known equations in suitable limits: the quantum linear Boltzmann equation if the internal degrees of freedom are 
neglected |3] , the master equation for an immobile test particle if the translational degrees of freedom are not relevant 
I14j , as well as quantum classical Markovian master equations if one or both kind of degrees of freedom can be 
treated as classical. Note that the natural bases in the derivation was given by momentum for the motional degrees of 
freedom and energy for the internal ones, the latter corresponding to the channel basis of scattering theory. In these 
cases different channels are only coupled through the collision term. If another internal basis can be of interest, also 
coherent tunnelling effects appear. 

We have further focused on the situation in which the internal degrees of freedom, in spite of influencing the collisional 
scattering cross section, are not probed by the measuring apparatus and therefore have to be averaged out from the 
set of the observed dynamical variables, thus effectively becoming part of the environment. The equation obtained in 
this situation is no more of Lindblad type, but rather takes the form of a generalized Lindblad structure [HIE]- It can 
therefore describe behaviour quite different from that characterizing a Markovian dynamics. Solving these equations 
in the position representation for an initial factorized state, we have obtained an explicit expression for the visibility 
reduction in interferometric experiments when internal degrees of freedom are involved. The behaviour of the visibility 
can indeed be quite different from the exponential decay corresponding to a Markovian dynamics, showing up e.g. 
oscillations and revivals. 

The interplay between different degrees of freedom in a bipartite system is a natural source of non-Markovian be- 
haviour, when either degrees of freedom cannot be controlled, thus acting as an environment. This scenario has here 
been studied in a concrete setting, assuming a factorized initial state, and describing the dynamics in terms of a gen- 
eralized Lindblad structure. Such a choice of initial condition is relevant for the considered interferometric setting, it 
would be however of great importance to consider initially correlated states, which naturally appear when considering 
a non-Markovian dynamics in a strong coupling regime. It would further be of interest to study whether, instead 
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of using this generalized Lindblad structure, one can obtain a closed description of the reduced system dynamics in 
terms of a master equation with a memory kernel, at least in some simplified situations. We plan to address these 
topics in future research work. 
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Appendix A: Microscopic derivation of the master equation 

We here address the derivation of the Markovian master equation for the description of the dynamics of the test 
particle starting from the general expressions Eq. ([s]) and Eq. Q. The scattering and rate operators appearing in 
these equations are best expressed using the factorization of the total Hilbert space "Htot = 'Hsys <8) "Hgas ='Hcm'8)'Hrei, 
according to 

|P,j)(P',j|®|p)(p'|gas = |rel(p,P),*)(rel(p',P'),j|rci®|P + p)(P'+p'|cm, (Al) 

where the Hilbert space (D" associated to the internal degrees of freedom is part of Hre\ [22], the notation is the 
same as in Sec. |ll] In fact both these operators act in a trivial way on centre of mass coordinates: f = Icm ^ Tq and 
T = Icm i8) Tq. The operator Fq is given by 



= !^^|rel(p,P)|atot(rel(p,P),j)®|j)(j|, (A2) 



where (Ttot(rel(p, P ),_;/) is the total cross section, depending on initial relative momentum and internal state. The 
relation W\ 



P.) (A3) 



links the operator Tq to the multichannel complex scattering amplitudes fkj{Pf iPi)^ referring to scattering from an 
initial momentum pi and internal state j to a final state with momentum pf and internal state k. According to 
standard usage in scattering theory we call channels the asymptotically free internal energy eigenstates of the system. 
The differential cross section is given by akj{pf,Pi) = (|p/|/|pi|)|/fcj(p/,Pi)P, so that the total cross section appearing 



in Eq. (A2) reads (Jtot{Pi,j) ^J2kl dPf^kj{Pf,Pi)- 



1. Evaluation of the Lindblad structure in momentum and internal state basis 

a. Incoherent contribution 

We now first concentrate on the evaluation of the contribution given by Eq. ([s]), which under suitable approximations 
can be cast in Lindblad form, closely following [T^ , where the special case of a test particle without internal structure 
was dealt with. To this end we consider the matrix elements of Cp in the momentum and channel basis {|P,«)} 
of the Hilbert space L^(R'^) (g) (D" associated to the test particle. Denoting by /i(p) the stationary gas momentum 



distribution and exploiting the relations Eq. ( A2 ) and Eq. ( A3 ) for the relevant operators we can express the result 
as is Eq. ([s]), where the complex rate functions 

Mll{P,P',Q) j dpoA*(Po)(rel(po-Q,P),*|Tory'|rel(po,i'-Q),j) 

x(rel(po,P' - Q) iry'TjreKpo - Q,P')) (A4) 

have been introduced, and |ri| denotes the volume in which the gas is confined. Note that M|/(P, P;<5) can be 
interpreted as classical rates for scattering of the test particle with momentum P ~ Q and internal energy eigenstate 
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j to a final state with momentum P and internal energy eigenstate i. The contributions at r.h.s. of Eq. ^ therefore 
play the role of quantmxi gain and loss term also depending on the internal degrees of freedom involved. 

Relying on Eq. (5| we can now deal just with the complex rates M^I,{P, P'; Q) defined in Eq. ( A4|. To proceed, it 
is helpful to introduce the following functions of po 



Pi 



rel 



P + P' 



POi 



Pf = rcl [po-Q, 



p' 



(A5) 



which denote the mean of the pairs of initial and final relative momenta appearing in M^(P,P';Q) and which are 
related by pi — Pf = Q. Introducing also q — rel (0, (P — P') /2) the complex functions M^^j^{P , P' \Q) can be 
expressed as an average over the gas distribution function /i of a complex density in the centre of mass frame 



Mil -P'; Q) = / dpoKPo)mli {pf,p^; q 



with 



Jl ( \ (27^^)% , •i-r ri/2| , A 

mik{pf,Pi;q) = —^{Pf + Q,A^o^Q \Pt + q,j) 

-l/2Tt 



(A6) 



Pt - q,i\^o '^llPf - q,k). 



(A7) 



Evaluating this formula with the expressions Eq. ( A2 ) and Eq. (LA3|) for Fq and Tp respectively, we obtain 



TO-i {Pf,P^;q) 



^^rl^\p. + q,jK^'ip.-q.f) 



2m* 



xS 



' Pi , [pf -Pi) q 



2m. 



Si 



'xfijiPf + qlP^ + q)fki{pf q^Pi q), 

i-gas IpI CTtot (p , j) I Tdt, is the eigenvalue of the operator Fq relative to the state |p, j) 



(A8) 



where ro(pi, j) = 

The expression given by Eq. ( A8 ) looses its meaning in the infinite volume limit, due to appearance of the arbitrarily 



large normalization volume. This point has been extensively discussed in [12 . It is to be traced back to the fact that 
the operator F in order to provide the actual rate of collisions should involve a projection on the subspace of incoming 
wave packets, which is not accounted for in Eq. (A7). To do this, we are going now to evaluate the operator F on a 
properly modified state of the relative motion. 



Before that, it is convenient to focus our attention on the two delta- functions appearing in Eq. (A8|: employing 
the relation (5(a)(5(6) = 2^(a + 6)<5(a — &), we can rewrite them as the product 



1^ / P)-P\ + Eu 

2 I 2m* 2 



(Pf-Pi ) ■ q . 



kl 



These two constraints ensure that the scattering amplitudes appearing in Eq. (AS I are evaluated on shell. The 



function mff^ {Pf,Pi\q) gives a significant contribution to the integral in Eq. ( A6) when the two energy differences are 



approximately equal, so that £ij = £ki, leading otherwise to rapidly oscillating phases, and this is actually a necessary 
condition in order to obtain a completely positive time evolution 13 . This implies in particular that integrating the 
generalized function m^|. {pf,Pi; q) with a function g{q), the contributions deriving from the parallel component of q 
vanish 



k iPf,Pi.q)9{q) 



{Pf,Pt;q±)giq±) 



(A9) 



We now therefore evaluate m^|. {Pf,Pi\qA^) with a properly modified state of relative motion, which takes into account 
the restriction of the expression to states which actually describe a colliding pair. To this end we write the complex 



rate m^], as 



(p/ + q_L,i|ToFy^exp (^i — 



pp. ® \3){l\ 



X exp 1 1 



X'ci • g-LA rl/2Tt 
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where pp. denotes an improper state of relative motion 



Pp. 



Pi){pi 



(AlO) 



Since the rate operator V should have a vanishing expectation value for those states of the relative motion that are 
not of the incoming type, we make the replacement 



PP.®\3){1\^ P'p^®\j){l\ = 



dxn 



dx 



-Lp. 



Pt 



w 



3){Pr 



W 



,l\. 



(All) 



This corresponds to a restriction of the Wigner function associated to the improper state of relative motion Eq. (AlO) 
from the entire normalization volume to a cylinder pointing in the direction pi, with base surface Sp^ and height 
Ap; . As in the case without internal degrees of freedom Ap. is approximately the distance travelled by the particle 
between two subsequent collisions, while Sp^ can now be taken as the geometric mean of the total cross-section of 
the involved channels [5], that is to say |Sp. | = \/o'(pi, j)cr(pj, /) and |Ap. | = {\pi\At)/m^,, with At the typical time 
interval between two subsequent collisions. 

Putting the new state Eq. ( All ) into the equation Eq. ( A7) and using the expressions of the matrix elements of Tg 
and Fq, we get 



dxn 



|ApJ 



p)-pi 

2m* 



x6 




{Eki — 



Pt + q± 



w 

~2 



Pt (1± 



w 

~2 



w 



^\ <7 [Pi + q± 



Pi (1± 



w 

~2 



where we have exploited once again the relation 6{a)d{b) = 2S{a + b)6{a — b). We now first perform the integral over 
t(j||p. , which denotes the component of w parallel to Pi, thus evaluating the second delta-function at r.h.s. of the 
previous equation, so that the dependence on x^^p. only appears in the term 

dxwr,. 



|Ap. 



\Pt 



-{8ij -8ki)] p 



The phase of the integrand varies very quickly for {£ij — £ki) ^ h/At, where At is the typical time elapsing be- 
tween collisions, so that as already discussed its contribution vanishes unless £ij — Ski, corresponding to a rotating 
wave approximation, assuming a separation of time scales between internal and translational dynamics [3]. Further 
considering the integral over x±p. as an approximate expression for 6 {w±p. ) we are led to 



iPf,Pi;Q±) 



^gas jl 



XikfijiPf + Q^^Pt + Qi-) 



'Jki {Pf - g-L,p« -Qi-)s{ + £, 



\Pi 



v/f7(Pj,j)cr(p»,0 



where the xH ^'^^ ^^^^ ^ Kronecker's delta factor, being defined according to 

1 if — £ki 
otherwise 
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In the last two terms we can disregard the dependence on gj_ because we expect that a qj_-integration will average 
out the "far off-diagonal" contributions with large modulus |g_L|, where the phases of the two scattering amplitudes 
are no longer synchronous. In conclusion, we have 



II 2 



Ei — Ei 



'Jij {Pf + q±,Pt + Q±) fki (Pf - Q±,Pi - 



(A12) 



This relation determines the complex rate functions {P,P';Q) through Eq. (A6|, and therefore the dissipative 
part of the master equation according to Eq. ([s]). 



b. Energy shift 



As a last step in the determination of the structure of the master equation we need to evaluate the contribution 
given by TZp. In the same notation as above, and within the same approximations, we directly obtain 



{P.,z\np\P\k) = {P,t\t Trgas f fri/2Re(T)ri/2, 



lp',fc) 



1^ (£;;■'■ (P)(p,j|p|p',fc) - Ei'^{p'){p,^p\p',j)) , 



(A13) 



with 



E^{P) = -27Tte'^xit I dpoKPo)Re (rel (po, P) , rcl {po,P))] • 
It is worth noting that for the case of a non-degenerate free internal Hamiltonian this formula reduces to 



(A14) 



1 



P,i\lZp\P\k) = - {EliP) ^ E^^iP')) {PAp\P\k) 



with 



K{P) 



rfPoM(Po)Re [J,i (rel (po, P) , rel (po, P))] ■ 



(A15) 



2. Operator expression of the master equation 



We now recast the master equation Eq. (jlj, whose matrix elements are given by Eq. ^ and Eq. {H^ , m a way 
which allows to express it in a representation- independent form. The key point is to show that M^j, {P,P';Q) ~ 
J dpop{po)mll. {pf,Pi; q) can be factorized into two terms, one depending on P and the other on P' . 

Changing the integration variable from po to pi and using the relations Eq. ( A5 ) to obtain po = Pi + {Pf + P) rri /M+ 
qm/msf = Pi + {pf + P')m/M — qmjm^^ we have 



,1/2 



P. + ^{P!^P) 



-qi- 



-(p,+P')--.. 



'^hj {Pf + Q±,Pt + fki (Pf - Q^-^Pi - : 



where we replaced q by qj^ in the arguments of /i^^^, in accordance with Eq. ( A9|. Remembering that Pi ~ Pf — Q 
and q = to* {P — P') / (2M), we consider the change of variable 



m 

Pi — Pi 

TO* 



TO P_L + P| 

M 2 



f. 



TO, Q 

TO* 2 Q^/m 



Q 



(A16) 



to obtain the desired factorization. If we further consider that the delta function 6 (p • Q/m) restricts the p- integration 
to the plane Qj^ = {p G K'^ : p • Q = 0} we finally arrive at the expression Eq. ([6]), with Lij (p, P; Q) as in Eq. 
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which allows to obtain the operator expression of the master equation given by Eq. (|9|) and Eq. ( |10[ ). If the gas 
distribution function fi is given by a Maxwell-Boltzmann probability density ^^(p) — l/ijT'^^^p^-j) exp(— p^/p|), where 

Pf} = y^2m/ P is the most probable momentum at temperature T — l/(fcs/3), these functions can be expressed in 
terms of the dynamic structure factor for a Maxwell-Boltzmann gas [28l [29] 



In fact, using the relation 




SMBiQ,E) = d^-ex^pl--^ '^ '-]. (A17) 



qM. [p. + + (1 + f + g^Qj = M,(P.)^MB(Q,i?(Q,P) +^.,), 

with E{Q, P) =: {P + Q)^/2A/ — P^/2A/ the energy transferred to the centre of mass in a collision changing the 
momentum of the test particle from P to P + Q, Eq. Q can be written as [6] 



L,, (p,P;Q) - J'^Mp±)\/Smb{Q,E{Q,P)+£,,) 

X/. (rel (p„ P,) - f + rel (p„ P,) + | + 

(A18) 

In the latter expression for the Lindblad operators, the dynamic structure factor appears evaluated for an energy 
transfer corresponding to the sum of the contributions for centre of mass and internal state, as naturally expected. As 
discussed in [4j, the dynamic structure factor describes momentum and energy transferred to the test particle when 
scattering off a macroscopic system, thus allowing for a more transparent physical understanding of the structure of 
the Lindblad operators. 
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